Experimental signature of phonon-mediated spin relaxation 
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We observe an experimental signature of the role of the phonons in spin relaxation between triplet 
and singlet states in a two-electron quantum dot. Using both the external magnetic field and the 
electrostatic confinement potential, we change the singlet-triplet energy splitting from 1.3 meV to 
zero and observe that the spin relaxation time depends non-monotonously on the energy splitting. 
A simple theoretical model is derived to capture the underlying physical mechanism. The present 
experiment confirms that spin-flip energy is dissipated in the phonon bath. 

PACS numbers: 03.65.W, 03.67.Mn, 42.50.Dv 



Relaxation properties of a quantum system are 
strongly affected by the reservoir where energy is dis- 
sipated. This has been seen clearly for electron spins 
embedded in nanostructures. Spin relaxation times T\ 
up to a few ns have been observed for free electrons in 
a two dimensional electron gas (2DEG) where energy is 
easily given to the motion In quantum dots, the 
discrete orbital energy level spectrum imposes other en- 
ergy transfer mechanisms. Near zero magnetic field, the 
electron spin can directly flip-flop with the surrounding 
nuclear spins, inducing short Ti's of the order of /xs 0- 
When a small magnetic field is applied, direct spin ex- 
change with nuclei is suppressed. Lattice vibrations, i.e. 
phonons, are expected to become the dominant reservoir 
in which spin-flip energy can be dissipated. This dissipa- 
tive mechanism is inefficient and very long spin relaxation 
times follow HHSH0. 

For spin relaxation involving phonons, two physical 
processes are important. The coupling between electron 
orbitals and phonons is responsible for dissipation of en- 
ergy and the spin-orbit interaction provides the essential 
mixing between different spin states so the spin states 
are coupled through electron-phonon interaction 
Energy conservation requires that the phonon energy cor- 
responds to the energy separation between the excited 
and the ground spin state. Changing the energy separa- 
tion affects the efficiency of the electron spin relaxation 
in two ways. First, since the phonon density of states in- 
creases with energy, the relaxation rate is expected to 
increase with energy as well. Furthermore, since the 
electron-phonon interaction is highly dependent on the 
matching between the size of the dot and the phonon 
wavelength [lfilllj . we expect a suppression of relaxation 
for very large and for very small phonon wavelengths 
in comparison to the dot size. Mapping the relaxation 
time T\ as a function of the energy splitting between the 
two spin states will provide insight in both the electron- 
phonon interaction and the spin-orbit coupling as well 
as an understanding of the limitations on T\ . This is of 
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FIG. 1: a) Dependence of the energy splitting AEs.t on to- 
tal magnetic field B. (b) Dependence of the energy split- 
ting AEs.t on the voltage Vr applied on gate 'T' at B= 0. 
Inset: Scanning electronic micrograph showing the sample 
design. The 2DEG, located 90 nm below the surface of a 
GaAs/AlGaAs heterostructure, has an electron density of 
1.3 ■ 10 15 m~ 2 . By applying negative voltages to gates L, 
M, T and Q we define a quantum dot (white dotted circle) 
and a QPC. Gate P is used to apply fast voltage pulses that 
rapidly change the electrochemical potentials of the dot. We 
tune the dot to the few-electron regime |T^|. and completely 
pinch off the tunnel barrier between gates L and T, so that 
it is only coupled to one electron reservoir at a time (lt^ . A 
voltage bias of 0.7 mV induces a current through the QPC, 
Iqpc, of about 30 nA. Tunneling of an electron on or off the 
dot gives steps in I QPC of 300 pA [H[]]l. The QPC mea- 
surement bandwidth is 100 kHz. 

particular relevance in the context of both spintronic and 
spin-based quantum information processing devices [l2| . 

Here, we study the spin relaxation time from triplet 
to singlet states for different energy separations in a sin- 
gle quantum dot containing two electrons. Singlet and 
triplet states have respectively two electrons in the low- 
est orbital and one electron each in the lowest and in the 
first excited orbital. In the experiment, the energy split- 
ting AEs,t between these two-electron spin states could 
be tuned from 0.9 meV to zero with a perpendicular mag- 
netic field and from 0.9 meV to 1.3 meV by deforming 
the dot potential [Hlli|. 

All the experiments are performed in a dilution refrig- 
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FIG. 2: (a) Voltage pulses applied to gate 1 P' for the relax- 
ation measurement. The starting point is a dot with one elec- 
tron in the ground state (initialization). During the pulse, 
the singlet and triplet electrochemical potentials are below 
the Fermi energy and a second electron tunnels into the dot. 
Due to the difference in tunnel rates 6], most likely a triplet 
state will be formed. We allow relaxation to occur during a 
waiting time that we vary. After the pulse, both electrochem- 
ical potentials are moved back above the Fermi energy and 
an electron tunnels out. This last step allows us to read-out 
the spin state, (b) Schematic of the AIqpc induced by the 
voltage pulse on gate ' P' . If the state was singlet, a step 
from a slow tunneling event is added to the QPC response 
just after the read-out pulse. If the state was triplet, the tun- 
neling event is too fast to be observed, (c) After averaging 
over many single traces, a dip is observed and its amplitude 
is proportional to the probability of having singlet present in 
the dot. (d) Relaxation curve obtained for B — 1.02 T con- 
structed by plotting the dip amplitude of the averaged traces 
at a pre-defined time after the read-out pulse. The relaxation 
time, Ti = 0.79 ± 0.05 ms, is extracted from an exponential 
fit to the data (all the data are taken with a 100 kHz low- 
pass filter). Inset : curve resulting from the averaging over 
500 individual traces for the longest waiting time (20 ms) and 
for the shortest waiting time (300 /is), offset by 100 fis and 
0.2 nA for clarity. 



erator with a quantum dot and a quantum point contact 
(QPC) denned in a 2DEG (see inset of Fig. []]»). The 
conductance of the QPC is tuned to about e 2 /h, making 
it very sensitive to the charge on the dot 01 • The sample 
is mounted at an angle cj) = 68° ± 5° with respect to the 
direction of the magnetic field B where ef> is derived from 
Shubnikov-deHaas oscillations. The magnetic field com- 
ponent perpendicular to the 2DEG is equal to B cos cf> 
(~ 0.385). The electron temperature was measured to 
be 180 mK from the width of the Coulomb peaks. The 
lattice temperature was 50 mK. 

We extract experimentally the energy splitting AEs.t 
between the singlet and the triplet states as a function 
of both magnetic field B and the confinement potential 
using a pulse spectroscopy technique 0. The depen- 
dence of AEs,t ° n B is presented in Fig. ^a). Up to 
0.4 T, AEs.t does not vary significantly with magnetic 
field which we relate to the elliptic nature of the dot at 
zero magnetic field [2jj. For B larger than 0.4 T, AE$ t T 
decreases, to a good approximation, linearly with mag- 



netic field. For energy separations below 100 fieV, the 
thermal broadening of the reservoir prevents us to mea- 
sure AEs,t- From extrapolation of the data, we can de- 
termine the magnetic field needed for singlet and triplet 
energy levels to cross: 2.82 ± 0.07 T. 

We measure the relaxation time for varying AEs.t- To 
be able to measure T\ close to the degeneracy point, we 
use a tunnel-rate selective read-out procedure (TRRO) 
@ (see Fig.|2). The measured spin relaxation time T\ as 
a function of B is presented in Fig. [3] The shape of the 
Ti dependence on magnetic field exhibits a striking non- 
monotonous behavior. From 0.4 T to ~ 2 T, correspond- 
ing to a decrease in the energy splitting from 0.8 meV to 
0.2 meV, the relaxation time first decreases, reaching a 
minimum of 180 /is. In between 2 T and the degeneracy 
point (2.82 T), T\ increases whereas the energy splitting 
continues to decrease. 

As a complementary study, we change AEs.t hi a dif- 
ferent way by controlling the electrostatic potential of the 
dot via the voltage Vr applied to gate 'T' and again look 
at Ti. The dependence of AEs.t on Vr is presented in 
Fig. IHb). With this second experimental knob, AEs.t 
can be varied from 0.8 meV to 1.3 meV. We interpret 
the change in the observed energy splitting as a conse- 
quence of a change in the dot cllipticity. A more positive 
Vr implies a more circular dot and a larger energy split- 
ting. We observe that T\ further increases with AEs.t 
as Vt is varied at B = T (see the inset of Fig.[3J). The 
maximum energy splitting reached at -530 mV, 1.3 meV, 
corresponds to a maximum of T\ = 2.3 ms. With both 
experimental knobs, we observe that when AEs.t is con- 
stant, T\ is constant too (respectively for Vt < —650 mV 
and B < 0.4 T). These observations clearly indicate that 
the most important parameter for the variation in the 
triplet-singlet relaxation time is their energy separation. 

The observed minimum in T± is precisely what one 
would expect for energy relaxation mediated by the 
electron-phonon interaction 0, Indeed, the en- 

ergy splitting AEs t T determines the relevant acoustic 
phonon energy (acoustic phonons are the only avail- 
able phonons for the explored energy range). At B 
~ 2 T, AEs.t ~ 0.3 meV, the associated half-wavelength, 
approximately 30 nm (the group velocity for acoustic 
phonons c s ~ 4000 m/s), is comparable to the expected 
size of the dot and therefore the coupling of the electrons 
in the dot to phonons is strongest. For energy separations 
smaller (larger) than 0.3 meV, the phonon wavelength is 
larger (smaller) than the size of the dot, the coupling 
to the orbitals becomes smaller and T\ increases. Taken 
together, all these observations strongly suggest that the 
phonon bath is the dominant reservoir for dissipating the 
spin-flip energy during relaxation. 

In order to get more insight in the role of the phonon 
wavelength, we present a simplified model of the energy 
relaxation process between triplet and singlet as a func- 
tion of their energy splitting AEs.t- From Fermi's golden 
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constructed from the lowest energy orbital and the first 
excited orbital (even though the contributions to triplet- 
singlet relaxation from higher orbitals can in fact be im- 
portant HH). In the notation |T_), |T+), |T ) and \S), 
both the orbital part (assuming Fock-Darwin states) and 
the spin part are present. Finally, we also neglect higher 
order (e.g. two-phonon) processes, which are important 
at small magnetic field |22|. 

In contrast to the one electron case 0,[_]], the spin-orbit 
interaction admixes directly the first excited states \T±) 
with the ground state \S). Due to the selection rules of 
the spin-orbit interaction, it does not affect \Tq) in lowest 
order j^. As a consequence, the spin relaxation time of 
\T ) can be much longer than \T±) 24]. However, we 
do not observe any signature of a slowly relaxing compo- 
nent in the experiment |25j . Since the spin-orbit coupling 
strength Mso is small in comparison with AE$\t (in the 
range accessed in the experiment), we can approximate 
the new eigenstates of the system as: 




FIG. 3: The spin relaxation time T_ as a function of the total 
magnetic field. The magnetic field where singlet and triplet 
states are degenerate is indicated by the dashed line. A mini- 
mum in Ti is observed around 2.2 T. The error bars represent 
70% confidence intervals. For energy separations close to de- 
generacy, the sensitivity of the measurement is reduced and 
the uncertainty in T_ increases. Inset: dependence of the re- 
laxation time Ti on Vr at B = T. 



rule, the relaxation rate between the triplet and the sin- 
glet states with energy separation AE$,t is proportional 
to their coupling strength through electron-phonon inter- 
action and to the phonon density of states at the energy 
AEs.t 11]- To obtain a simple analytical expression, we 
assume that the only effect of the perpendicular magnetic 
field, the Coulomb interaction between electrons and the 
modification of the potential landscape is to change the 
energy splitting. Especially, their effects on the spatial 
distribution of the wavefunctions are neglected and we 
neglect the Zeeman energy. Furthermore, we restrict the 
state space of the analysis to |T_), \T + ), \T ) and \S) where H, 
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FIG. 4: Relaxation rate as a function of the energy splitting 
AEs,t deduced from the experimental data. The circles and 
the triangles correspond to the experiment where we vary 
respectively the magnetic field and the dot potential. The 
solid (dotted) line is the curve for Mso = 0.37 ^teV (Mso = 
2.31 AteV) obtained from the simplified model. 
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In general, Mso is dependent on the magnetic field, 
but to simplify the discussion, we neglect this dependence 
[2ll |22| . Since the electron-phonon interaction preserves 
the spin, the coupling between \T±) and \S') has the fol- 
lowing form: 



{T' ± \H e , p \S') = 



so 



S,T 



((S\H e JS) - (T±\H e jT ± )) 



giq.n_j_giq.r2 j g j-j^g interaction Hamiltonian 
between electrons and phonons, q the phonon wavevec- 
tor and ri the positions of the electrons. One can then 
interpret the coupling between \T±) and \S') as the differ- 
ence of the electron-phonon interaction strength for the 
corresponding unperturbed states \T±) and \S). If the 
phonon wavelength is much larger than the dot size, the 
coupling to the phonons is the same for both states and 
the two terms will cancel. If the phonon wavelength is 
much shorter than the dot size, the coupling is small for 
each state separately. 

To provide a quantitative comparison to the data, we 
need to model the electron-phonon interaction. Following 
IB EH- w e assume bulk-like 3D phonons. For the energy 
separations discussed in our experiment, only acoustic 
phonons are relevant. The Hamiltonian H ep has then 
the following expression: 
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F z {qz) 
y/2pqc~Jh 
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where (q, j) denotes an acoustic phonon with wave vec- 
tor q = (qiucfe), j the phonon branch index and p — 
5300 kg/m 3 is the density of lattice atoms. The fac- 
tor F z (q z ) depends on the quantum well geometry and 
is assumed to be 1 in our model The speed of 

sound for longitudinal and transverse phonons are re- 
spectively ci — 4730 m/s and c t = 3350 m/s. We 
consider both piezo-electric and deformation potential 
types of electron-phonon interaction. In the consid- 
ered crystal, the deformation potential interaction is 
non-zero only for longitudinal phonons (with a coupling 
strength H = 6.7 eV). In contrast, all phonon polariza- 
tions j are important for piezo-electric coupling. The 
coupling strength depends on 9, defined as the angle be- 
tween the wavevector and the growth axis, and varies 
for different polarizations as e/3 Jiq = Aj{6)ej3 where 
ej3 = 1.4 x 10 9 eV/m @, HI- Due to the different de- 
pendence on q for both mechanisms (^fq for deformation 
potential interaction, l/^/q for piezo-electric interaction) , 
the piezo-electric (the deformation potential) coupling 
between electrons and phonons is dominant for energy 
separations below (above) 0.6 meV. From direct applica- 
tion of Fermi's golden rule, we derive the following ana- 
lytical expression for the spin relaxation rate 1/71 : 

AE| t q 2 sin 2 e 

1 /rp M| a 4 /H 2 A_E| T rir/2 , q . 5 „ ^272 

1 / T i=a£far{ h*cf Jo d6sm9e ^ 

+ Ej rj 2 de\Me)\ 2 Bin 5 e e ^ ) 

where a is the dot radius (in our model a is independent 
of A.Est an d is estimated to be 23nm, from the mea- 
sured single particle level spacing). This simple model 
reproduces the most important feature in the measure- 
ments, which is that the coupling to the phonons vanishes 
for large and small energy separations and is strongest 
when the phonon wavelength matches the dot size (see 
Fig.gJ. 

The spin-orbit strength Mso appears in the expres- 
sion of 1/Ti only as a scaling factor. With a value 
Mso — 0.4 ^eV (corresponding to a spin-orbit length 
equal to %/ 2am* Mso ~ 50 A* 111 ); the model reproduces 
the peak amplitude of the data quite well (Fig. 0] solid 
line). However, this value for Mso is about six times 
smaller than the values reported in (the dotted 

line in Fig. 0] corresponds to the relaxation rate using 
this value of Mso in the model) . The discrepancy could 
be the result of the exclusion of higher orbitals and the 
magnetic field dependence of Mso in our model [2lL Eij ] . 
Again, we emphasize that both curves have a maximum 
corresponding to a phonon wavelength matching the dot 
size. 

For single electron spin states, comparable variations 
of T\ with the energy splitting are expected although 
direct spin-orbit coupling between Zeeman sublevels of 
the same orbital is zero. To maximize the relaxation 



time of electron spin qubits, one needs then to choose 
an energy separation between the spin states such that 
the corresponding phonon wavelength is different from 
the dot size. To complete our study of spin relaxation, it 
will be interesting to rotate the sample with respect to 
the magnetic field since the spin-orbit coupling strength 
depends on the angle between the crystallographic axis 
and the magnetic field 0, E<| . 
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